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I.  Introduction 

Reinforcement  of  polymers  using  organic  and  inorganic  particles  (fillers)  has  become  extremely  common  in  a 
variety  of  practically  important  applications.  Traditional  applications  in  this  context  belonged  to  the  “colloid”  or 
the  “composite”  regime,  where  the  filler  size  was  typically  larger  than  the  size  of  the  polymer.  Not  surprisingly, 
the  composite  limit  has  also  had  many  advances  in  theoretical  models  and  simulation  approaches  for  predicting 
the  equilibrium  and  dynamical  properties  of  such  mixtures.!-  However,  more  recent  developments  in 
nanotechnology  applications  have  involved  polymer- filler  mixtures  in  the  “nanoparticle  limit”  where  the  size  of 
the  polymer  is  comparable  to  or  larger  than  one  or  more  dimensions  of  the  filler.!— b!  In  such  materials,  termed 
polymer  nanocomposites  (PNCs),  uniform  dispersion  of  the  filler  particles  results  in  significant  interfacial 
contact  between  the  polymer  and  the  filler,  which  in  many  cases  leads  to  new  and  novel  properties  arising  from 
the  unique  synergism  between  materials. 

Many  researchers  have  demonstrated  the  potential  of  PNCs  for  a  variety  of  applications4510-1 1  For 
example,  incorporating  nanoscale  dispersions  of  layered  clay  platelets,  carbon  nanotubes  and  nanosized  silica 
particles  into  polymers  has  been  shown  to  enhance  both  the  amorphous  and  the  rubbery  modulii  of  the  bare 
polymer  matrix  by  as  much  as  an  order  of  magnitude.2’3’12"14  Gas  barrier  properties  of  butyl  rubber  latexes 
was  shown  to  be  reduced  by  almost  2  orders  of  magnitude  upon  incorporating  just  10  wt%  of  vermiculite 
fillers.—  Addition  of  only  5%  by  weight  of  clay  platelets  were  shown  to  reduce  the  fire  hazard  of  nylon-6 
polymer  by  around  60%.!C!Z  Similarly,  polydimethylsiloxane  (PDMS)  filled  with  around  10%  by  weight  of 
clay  platelets  displayed  a  140  °C  higher  decomposition  temperature  compared  to  the  pristine  PDMS 
elastomer.—  Mixtures  of  polymers  with  carbon  nanotubes,  semiconducting  particles  and  magnetic  particles 
have  been  demonstrated  to  possess  novel  electrical,  optical  and  magnetic  properties.7’18"23 

Development  and  application  of  nanoscale  multicomponent  materials  such  as  PNCs  confronts  a  huge 
parameter  space  involving  an  interplay  of  constituent  selection,  lubrication,  processing  and  performance.!!-  In 
this  article,  we  specifically  focus  on  an  issue  which  has  commanded  significant  attention  in  PNCs,  viz., 
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dispersion  control  of  nanoparticles  in  polymer  matrices.  Most  combinations  ofpolymers  and  pristine 
nanofillers  tend  to  be  immiscible,  with  the  fillers  undergoing  aggregation  either  due  to  strong  Van  der  Waals 
interactions  between  themselves  or  due  to  polymer-mediated  interparticle  attractions5,6’3’9’24"26  (an  example 
illustrating  the  dispersion  issue  is  displayed  in  Fig.  1  — ).  Filler  aggregation  usually  has  a  catastrophic  effect  on 
the  properties  of  PNC  s,  since  many  characteristics  of  interest  in  PNCs  typically  depend  on  the  significant 
interfacial  area  afforded  at  the  nanoscale  —  a  feature  which  is  considerably  eroded  due  to  the  aggregation  of 
the  fillers.  Consequently,  a  significant  issue  confronting  the  development  of  PNCs  for  applications  is  to  achieve 
well- dispersed  filler  configurations  within  the  polymer  matrix. 


Fig.  1  (Adapted  from  ref  27  with  permission).  TEM  and 
optical  micrographs  of  poly(methyl  methacrylate)  (PMMA)-C6o 
PNCs  at  different  Cgo  wt%.  The  dark  features  are  C(,o 
agglomerates.  At  C^q  wt%  less  than  0.01,  C6o  units  are  seen  to 
be  well- dispersed  with  agglomerate  size  units  of  scale  around 
20  nm  in  diameter.  At  C6o  wt%  =  0.01  nanoscale  agglomerates 
are  seen  to  coexist  with  micron  sized  agglomerates.  For  C$o 
wt%  =  0.05,  the  C50  are  seen  to  agglomerate  into  bundles  of 
micron  sizes. 


A  simple  understanding  of  the  issues  involved  in  the  context  of  achieving  particle  miscibility  and 
dispersability  can  be  obtained  by  considering  the  surface  tension  of  a  spherical  particle  in  a  polymer  melt 
displayed  in  Fig.  2.  It  can  be  seen  that  for  cases  where  there  are  unfavorable  interactions  between  the  polymer 
melt  and  the  particle,  the  surface  tension  is  positive.  One  would  expect  that  uniform  dispersion  of  particles 
would  not  be  favored  in  such  situations  and  that  the  particles  would  prefer  to  either  aggregate  or  separate  out 
of  the  matrix.  In  contrast,  favorable  polymer-particle  interactions  lead  to  a  negative  surface  tension.  In  such 
cases,  intimate  mixing  between  the  polymer  and  particles  lead  to  a  lower  free  energy  and  therefore  might 
constitute  conditions  favoring  dispersability  of  the  particles. 


Fig.  2  Surface  tension  )'  (in  units  nondimensionalized  with 
kBT!a 2 ,  where  kB  denotes  the  Boltzmann  constant  and  a 
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represents  the  Kuhn  segment  length  of  the  polymer)  between  a 
particle  of  size  R  (in  units  nondimensionalized  by  the  polymer 
radius  of  gyration)  and  a  compressible  polymer  melt  of  flexible 
chains.  The  parameter  A  (termed  the  adsorption  length) 
quantifies  the  adsorption  interaction  between  the  polymer  melt 
and  the  particle.  Explicitly,  A/a  =  Xs~  Xsc,  where  Xs  denotes 
the  repulsion  energy  per  monomer  (in  kBT units)  between  the 
surface  and  the  polymer.  Xsc  denotes  the  critical  interaction 
value  between  the  polymer  surface  and  the  particle  at  which 
the  adsorbed  surface  excess  of  the  polymer  is  identically  zero. 

Positive  (negative)  values  of  A  correlate  to  repulsive 
(attractive)  interactions  between  the  particle  and  the  polymer. 

The  polymer  melt  is  modeled  using  a  compressible  polymer 
melt  model  described  in  ref  28  with  the  nondimensional 
compressibility  parameter  chosen  as  0. 1 .  The  interfacial  tension 
was  determined  using  the  formalism  described  in  ref  29. 

A  number  of  experimental  efforts  have  used  the  above  guiding  principles,  and  have  focused  on  controlling 
the  state  of  dispersion  of  nanofillers  by  a  combination  of  one  or  more  of  the  following  strategies:  (i)  choosing 
the  polymer- filler  combinations  appropriately  to  take  advantage  of  favorable  polymer-filler  interactions;  (ii)  by 
“functionalizing”  the  fillers  or  the  polymer  matrix  with  a  variety  of  anionic  or  cationic  oligomeric  surfactants  and 
grafted  polymers  to  exploit  favorable  interactions  between  the  functionalizing  group  and  the  polymer 
matrix;3’9’30,25’26’31  and  (iii)  using  external  fields  such  as  electric,  magnetic  or  flow  fields  to  disrupt  the 
equilibrium  (aggregated)  state  and  maintain  it  in  a  nonequilibrium  dispersed  and/or  aligned  configuratioa— 

While  there  have  been  a  number  of  experiments  reported  along  the  above  lines,  the  outstanding  issue  lacing 
the  development  of  theories,  models  and  computer  simulations  to  aid  the  design  of  such  strategies  is  the 
following:  “for  a  specified  combination  of  matrix  polymer(s),  filler(s)  and  the  functionalizing  moieties,  can  the 
expected  structure  and  properties  of  PNC  dispersions  be  predicted?”  While  the  guiding  principles  discussed 
above  have  typically  served  as  a  good  qualitative  rule  of  thumb  for  choosing  the  different  components,  such 
simplistic  principles  suffer  from  several  limitations.  Explicitly,  (i)  such  considerations  do  not  account  for  the 
particle  concentration  effects  as  might  be  embodied  within  a  detailed  “phase  diagram"  for  the  polymer-particle 
mixture;  (ii)  the  influence  of  size  and  chemical  characteristics  of  the  functionalizers  (either  small  molecule  or 
polymeric)  is  not  accounted  for;  (iii)  dispersion  of  particles  in  structured  matrices  such  as  the  self- assembled 
phase  of  a  block  copolymer,  cannot  be  gleaned  from  the  above  considerations. 

Not  surprisingly,  overcoming  the  above  limitations  have  constituted  the  focus  of  several  recent  theoretical 
investigations.  Various  methods,  including  computer  simulations,33-35  liquid  state  theories—  and  polymer  field- 
theoretic  approaches'^  have  been  used  to  address  several  aspects  of  the  above  issues.  In  our  own 
research,  we  have  adopted  the  use  of  coarse-grained  modeling  approaches  to  delineate  the  mechanistic 
features  underlying  the  structure  and  dynamical  properties  of  nanoparticle-polymer  mixtures.39-47  Such 
coarse-grained  strategies  typically  involve  the  use  of  simple  micromechanical  models  to  represent  the  different 
components.33 ,48~52  For  instance,  a  common  approach  is  to  model  the  polymer  as  a  connected  sequence  of 
segments,  where  each  segment  is  understood  to  represent  a  collection  of  molecules  or  atoms.—  Moreover,  the 
polymer  chain  is  typically  (albeit,  not  necessarily)  assumed  to  be  fully  flexible  and  behave  as  an  elastic  spring 
or  a  Gaussian  coil  The  particle  fillers  in  such  models  may  be  represented  as  either  hard  spherical  or 
anisotropic  objects.  Moreover,  the  different  physicochemical  interactions  are  represented  by  ignoring  the 
specific  chemical  identities  of  the  monomers  and  the  resulting  detailed  interaction  characteristics.  Instead 
simpler  model  interaction  potentials  are  used  to  characterize  the  interactions  between  segments  of  the  polymer 
and  the  particle  fillers.  Such  simplifying  assumptions  have  enabled  the  implementation  of  analytical  theories 
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and/or  long  time  and  length  scale  simulations  which  allow  one  to  discern  the  equilibrium39-4 1  -43-44-47  and 
nonequilibrium45  46  structural  characteristics  of  the  nanofiller  dispersion  in  polymer  matrices.  We  note  that  the 
main  utility  of  such  coarse-grained  approaches  lies  in  their  ability  to  distill  and  characterize  physical 
phenomena  of  interest  in  terms  of  a  few  macroscopic  parameters. 

In  this  article,  we  review  some  of  our  recent  contributions  using  such  coarse-grained  models  and 
simulations  to  address  the  physical  principles  underlying  the  equilibrium  aspects  of  dispersion  strategies. 
Specifically,  we  focus  on  three  issues:  (i)  dispersion  and  phase  behavior  of  nanoparticles  in  homopolymer 
matrices;39-41 -43-44-47  (ii)  dispersion  in  mixtures  of  homopolymers  with  grafted  nanoparticles;  (iii)  self  assembly 
and  organization  of  nanoparticles  in  block  copolymer  matrices.—  Each  of  these  topics  serve  to  illustrate  that 
the  dispersability  of  nanoparticles  may  exhibit  far  more  complexities  than  the  simple  miscibility  trends  one  may 
deduce  using  the  results  ofFig.  2.  Considering  the  breadth  of  this  field  and  the  rate  at  which  new 
developments  are  reported,  no  attempt  is  made  to  render  this  a  comprehensive  review  of  each  of  the  topics. 
Consequently,  our  review  is  very  selective  and  focuses  only  on  issues  of  our  expertise  and  specifically  on  the 
developments  which  have  arisen  out  of  our  own  research.  However,  where  appropriate,  related  theoretical 
and  experimental  studies  and  some  possible  future  directions  are  briefly  mentioned. 

II.  Phase  behavior  of  nanoparticles  in  homopolymer  matrices 

The  simplest  model  system  pertinent  to  nanocomposite  materials  is  that  of  a  mixture  of  a  homopolymer  matrix 
(melt  or  solution)  with  nanoparticle  fillers.—  Whence,  the  first  class  of  studies  we  review  pertain  to  our 
contributions  toward  the  modeling  of  the  dispersion  characteristics  of  spherical  nanoparticle  fillers  in 
homopolymer  solutions  and  melts.  In  this  regard,  we  undertook  several  studies  with  an  objective  to  obtain  a 
fundamental  understanding  of  the  manner  in  which  the  polymer- polymer,  polymer- filler  and  filler-filler 
interactions  control  the  phase  behavioral  characteristics  of  such  systems.39’41-43’44-47  Specifically,  our  research 
was  focused  on  understanding  the  influence  of  particle  curvature  and  the  specifics  of  the  “nanoparticle  limit” 
upon  the  overall  phase  behavior  and  particle  structure  in  such  polymer- nanoparticle  mixtures. 

That  interesting  curvature  effects  may  manifest  is  clearly  evident  even  from  Fig.  2  where  it  is  seen  from  the 
size  dependence  of  the  surface  tension  that  the  unfavorable  or  favorable  mixing  effects  tend  to  be  diminished 
for  smaller  sized  particles.  An  outstanding  question  then  is  “are  polymer-nanoparticle  mixtures  always  miscible 
or  is  there  a  potential  for  richer  phase  behavior  characteristics?  ”—  To  address  this  issue,  our  research  has 
used  coarse-grained  polymer  field  theories^^  to  address  the  dispersion  characteristics  of  nanoparticles 
mixed  with  polymer  solutions.  Our  approach  was  based  on  a  rigorous  formalism  which  integrates  out  the 
polymer  “degrees  of  freedom”  to  derive  polymer- mediated  effective  interaction  potentials  between  the 
nanoparticles.  Such  interaction  potentials  were  then  combined  with  the  bare  particle-particle  interactions 
within  thermodynamic  theories  and/or  computer  simulations  to  shed  light  on  the  dispersion  and  structural 
characteristics  of  the  nanoparticles  in  the  polymer  matrix.  In  the  following,  we  first  briefly  outline  the  theoretical 
formalism  we  used  and  then  subsequently  discuss  selected  results  from  a  few  of  our  applications. 

A.  A  mean-field  approach  to  the  phase  behavior  of  polymer-nanoparticle  mixtures 

We  developed  a  mean- field  theoretic  formalism  to  address  the  phase  behavior  of  polymer- nanoparticle 
mixtures.??^-!.  We  adopted  an  implicit  solvent  framework  comprised  of  a  two  component  system  of  just 
particles  (c)  and  polymers  (p),  interacting  with  each  other  by  effective,  solvent- averaged  interaction  potentials. 
Our  formalism  used  a  McMillan- Mayer  like  solution  theory  to  formally  recast  the  statistical  mechanics  of  such 
a  two-component  system  of  polymer  and  particles  into  a  single  component  system  of  just  particles  which 
interact  by  polymer- mediated  effective  interaction  potentials  in  addition  to  their  bare  interparticle  potentials. 
Explicitly,  in  a  grand- canonical  framework  for  the  polymer  solution,  the  polymer- mediated  pair- interaction 
potentials  between  nanoparticles  Z/(r„  rj)  can  be  shown  to  be  expressible  as/A/if. 
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U(rh  rj)  =  h S2(r/,  r/,zp,  V)  -  h Si(r (jzp,  V)-foZ 

i(tjiZp,  F)  +  h  £,o(zjp) 

In  the  above  equation,  £„(r„  i ),  . . r„;  zp)  in  general  denotes  the  grand  canonical  partition  function  for  the 
polymer  solution  at  a  fixed  activity  coefficient  zp  containing  n  particles  fixed  at  the  positions  r„  r;. .  ,r„. 

To  obtain  E„(r y,  i},  . . rm;  z/;)  for  the  polymer  solution,  we  implemented  a  mean- field  theoretic  approach 
commonly  known  as  polymer  self-consistent  field  theory  (SCFTf.EfM  In  a  nutshell,  polymer  SCFT 
enumerates  the  statistical  features  and  the  thermodynamics  of  an  interacting  system  of  polymer  chains  by 
mapping  them  onto  an  equivalent  system  of  noninteracting  chains  in  the  presence  of  pseudo  chemical  potential 
field(s),  which  in  turn  embody  the  interactions  of  a  specified  polymer  chain  with  other  polymer  chains.  The 
basis  for  such  a  formalism  is  established  through  field-theoretic  techniques,  which  can  be  used  to  demonstrate 
that  the  thermodynamics  of  the  system  of  noninteracting  chains  serves  as  a  mean- field  approximation  to  the 
thermodynamics  of  the  system  of  interacting  chains.—  In  a  mean- field  approximation,  the  intersegment 
interactions  themselves  are  a  function  of  the  inhomogeneous  densities  of  the  appropriate  components.  The 
latter  are  in  turn  themselves  influenced  by  the  chemical  potential  field(s)  acting  on  the  polymer  segments. 
Consequently,  implementation  of  SCFT  typically  requires  the  self-consistent  solution  of  a  set  of  field  equations 
for  the  chemical  potential  field(s).— 

To  implement  the  above  formalism,  we  adopted  a  commonly  used  model  for  the  polymers  termed  the 
Gaussian  thread  model.—  In  this  model,  the  monomeric  units  of  the  polymers  are  assumed  to  be  point-sized 
and  the  polymer  chains  themselves  to  be  elastic  threads  connecting  these  monomers.  The  polymer-polymer 
interactions  were  modeled  through  effective  excluded  volume  interactions  which  represent  the  combined 
effects  of  the  polymer-polymer  and  polymer-solvent  interactions.^!!  In  the  SCFT  formalism  for  polymer 
solutions,  the  effects  of  such  excluded  volume  interactions  are  replaced  by  a  spatially  inhomogeneous  chemical 
potential  field,  denoted  as  JT(r),  which  is  determined  as  the  solution  of: 

W(r)  =  BC4>  (r)  (2) 

where  B  represents  a  nondimensional  excluded  volume  parameter  and  C  denotes  the  nondimensionalized 
overall  polymer  solution  density.  The  field  4>(r)  represents  the  nondimensionalized  inhomogeneous  density  field 
of  the  polymer  segments.  The  above  equation  has  a  simple  physical  meaning  in  that  it  quantifies  in  a  mean- field 
manner  the  excluded  volume  interactions  experienced  by  a  segment  at  location  r  due  to  the  other  intra  and 
interchain  polymer  segments. 

Eqn  (2)  is  rendered  a  self-consistent  condition  by  requiring  that  the  density  field  <H,x)  be  itself  obtained  as  a 
result  of  the  statistics  of  the  noninteracting  chains  in  the  external  field  W(x).  The  main  feature  which  allows  for 
practical  implementation  of  this  formalism  is  the  feet  that  the  statistics  of  a  polymer  molecule  in  an  external  field 
can  be  determined  by  solving  a  diffusion  equation  for  its  distribution  functions.  For  instance,  for  a  polymer 
chain  under  the  action  of  an  external  field  W(r),  the  probability  q(  r,  s )  that  the  5th  segment  of  the  chain  lies  at 
the  position  r  (irrespective  of  the  starting  position  of  the  chain)  satisfies  :E!!! 

%  =  v*V q(r,s  =  0)  =  I  (3) 

Sr 

The  density  field  <;>  (r)  of  the  polymer  can  be  expressed  in  terms  of  q(  r,  5)  as: 

«*(r)  =  S\ds  q(r,s)q(r,N  -  s)  (4) 

where  N  denotes  the  number  of  segments  in  the  polymer.  Moreover,  the  grand  canonical  partition  function  of 
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the  polymer  solution  can  be  expressed  in  the  mean- field  approximation  as:- 

^  p/r  W1  +  7.  pr^(r, *  =  V )  (5) 

In  the  above  equation,  Z  represents  a  nondimensionalized  chemical  potential  for  the  polymer  solutioa 

Overall,  the  above  formalism  provides  a  means  to  obtain  the  grand  canonical  partition  function  Z.  In 
general,  the  above  formalism  also  allows  for  incorporating  a  variety  of  physical  polymer-particle  interactions. 
Typically,  such  interactions  manifest  as  appropriate  boundary  conditions  on  the  diffusion  eqn  (3).50’58,59  For 
instance,  a  purely  impermeable  wall  is  modeled  as  a  boundary  condition:  q( r,  s  >  0)  =  0  on  the  surface.  While 
effects  of  energetic  interactions  between  the  surface  and  the  polymer  can  be  incorporated  directly  as  an 
external  potential  on  the  polymer  segments,  situations  where  the  surface  exhibits  an  extremely  short  ranged 
interaction  (of  range  smaller  than  or  comparable  to  the  segment  sizes),  are  usually  modeled  through  a 
boundary  of  the  form: 


n-Vg(r,  5)  =  ~Aq 

on  the  surface.  In  the  above,  n  denotes  the  normal  to  the  surface,  and  A  1  is  a  positive  (negative)  length  scale 
quantifying  the  strength  of  attractive  (repulsive)  interactions.—  An  excellent  discussion  of  the  origin  of  boundary 
conditions  and  the  associated  numerical  details  can  be  found  in  the  monograph  by  Fredricksoa—  Solution  of 
the  SCFT  equations  with  such  boundary  conditions  allows  one  to  determine  Zn(r,-,  r y,  . . .,  r„;  zp),  which 
quantifies  the  grand  canonical  partition  function  of  the  polymer  solution  in  the  presence  of  fixed  particles.  In 
turn,  this  allows  one  to  deduce  the  polymer- mediated  interactions  through  eqn(l). 

In  general,  the  diffusion  eqn  (3)  does  not  admit  an  analytical  solutioa  Since  our  objective  in  this  research 
was  to  specifically  examine  the  phase  behavior  of  polymer-nanoparticle  mixtures  for  a  range  of  particle  sizes, 
especially  for  regimes  where  the  curvature  of  the  particle  proves  crucial  in  determining  the  physics,  we  solved 
the  diffusion  eqn  (3)  numerically.—  This  procedure  was  executed  in  spherical  coordinates  during  the 
computation  of  Si  (Z,  B )  and  in  bispherical  coordinates  during  the  computation  of  £2(1',,  r/,  Z,  B).  The  use  of  a 
bispherical  coordinate  system  allows  us  to  access  a  wide  range  of  particles  sizes  without  encountering  any 
artifacts  arising  from  geometrical  discretization  errors.  More  details  on  the  numerics  of  the  method  can  be 
found  in  Roan  and  Kawakatsu.—  In  the  next  section,  we  review  the  results  of  two  classes  of  studies  in  which 
we  have  used  the  above  formalism  to  model  polymer-nanoparticle  mixtures. 

B.  Nanoparticles  in  solutions  with  depleting  polymers 

Our  first  application  of  the  above  formalism  was  to  study  the  interactions  and  phase  behavior  of  nanoparticle- 
polymer  mixtures  for  the  case  where  the  polymers  and  particles  have  no  direct  energetic  interactions,  except 
insofar  as  the  polymers  being  only  excluded  from  the  interiors  of  the  particles  (commonly  known  as  the  “ 
depletion”  situation).—  We  note  that  this  model  is  the  simplest  among  the  class  of  models  characterizing  the 
behavior  of  polymer- nanoparticle  mixtures  in  which  there  are  unfavorable  enthalpic  interactions  between 
polymers  and  particles.—  In  such  cases,  it  can  be  expected  that  polymer  exclusion  from  the  particle  surfaces 
leads  to  an  effective  interparticle  attraction  between  the  particles  which  drives  the  aggregation  and  phase 
separation  of  particles.  Prior  theoretical  research^Ail  had  considered  the  magnitude  of  such  depletion 
interactions  in  the  “nanoparticle”  regime,  and  had  suggested  that  such  depletion  interactions  were  weak  for 
small  particles  and  that  nanoparticle-polymer  mixtures  may  be  expected  to  be  stable  against  demixing  arising 
from  such  depletion  attractions.  Our  research  was  specifically  motivated  to  analyze  these  curvature  effects  in 
more  detail  due  to  experimental  reports  which  contradicted  with  these  theoretical  considerations  and 
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demonstrated  that  in  the  nanoparticle  regime  smaller- sized  particles  may  actually  be  more  prone  to  phase 
separation  than  larger- sized  particles.— 

A  model  for  depleting  polymers  was  straightforwardly  implemented  in  the  context  of  the  formalism 
described  in  section  IIA  by  adopting  a  boundary  condition  q(  r,  s)  =  0  inside  the  surface  of  the  particles.—  This 
boundary  condition  embodies  the  inpenetrability  of  the  particle  surface  to  polymer  segments.  Using  such  a 
formalism,  in  ref  39  we  analyzed  the  phase  behavior  of  nanoparticle  dispersions  in  polymer  solutions.  A  first 
clue  towards  unraveling  the  peculiarities  of  the  nanoparticle  regime  arose  from  considering  the  thickness  of  the 
polymer  exclusion  (depletion)  zone  around  the  particles  (cf  Fig.  3a).ZUri  Explicitly,  the  latter  showed  that,  in 
the  nanoparticle  regime,  the  volume  of  the  polymer  depletion  layers  can  for  exceed  the  size  of  the  particles. 
Since  the  range  of  the  depletion  layers  are  also  expected  to  determine  the  range  of  the  polymer-mediated 
interparticle  potentials,  these  results  suggested  that  the  nanoparticle  regime  may  be  accompanied  by  significant 
multibody  interaction  effects.  Our  hypothesis  was  that  such  multibody  effects  may  render  the  effective 
interactions  and  phase  diagram  to  be  significantly  different  from  what  may  be  deduced  by  considerations  of  the 
infinite  dilution  situation  involving  a  single  particle. 


Fig.  3  (Results  adapted  with  permission  from  ref  39).  (a)  Depletion  thickness  of  polymer 
layers  (in  units  normalized  by  the  radius  of  the  particle  R)  as  a  function  of  the  radius  of  the 
particle  (in  units  normalized  by  the  unperturbed  radius  of  gyration  of  the  polymer  Rg).  4> 
represents  the  polymer  concentration  normalized  by  the  overlap  concentration,  (b)  Particle 
volume  fraction  dependencies  of  polymer- mediated  pair- interaction  potentials  U(r)  (ft  = 

\/kBT)  as  a  function  of  the  distance  between  the  particles  r.  The  polymer  concentration  was 
maintained  at  ^  =  0.2  and  the  particle  size  RJRg  =  5;  (c)  Fluid-fluid  coexistence  curves  for 
different  particle  sizes.  The  region  above  the  lines  represent  the  regions  of  immiscibility.  The 
dashed  lines  represent  the  tie- lines  for  RJRg  =  0.33. 

To  account  for  the  above  multibody  effects,  we  proposed  an  approximate  approach  within  the  pair 
interaction  formalism  by  rendering  the  chemical  potential  dependence  Z  of  the  polymer  concentration  C  to  also 
include  a  dependence  on  the  particle  volume  fraction.  The  latter  dependencies  were  in  turn  deduced  using  a 
free  volume  theory  model  for  binary  hard  spheres.—  The  influence  of  such  effects  are  displayed  in  Fig.  3b. 
which  displays  the  particle  volume  fraction  dependence  of  the  pair- interaction  potentials  (for  a  fixed  overall 
polymer  concentration).  It  is  clearly  seen  that  white  the  polymer-mediated  attractions  are  weak  for  dilute 
concentrations  of  particles,  the  polymer- mediated  attractions  become  much  stronger  for  nondilute 
concentration  of  particles.  We  translated  such  effects  into  a  phase- diagram  using  a  simple  thermodynamic 
perturbation  theory- like  approach  in  which  the  polymer- mediated  interaction  potentials  were  treated  as  a 
perturbation  to  the  bare  interparticle  hard  sphere  interactions.  From  the  results  displayed  in  Fig.  3c  we 
observe  that  nanoparticle-polymer  mixtures  do  indeed  show  a  large  region  of  immiscibility  over  a  wide  range 
of  particle- polymer  size  ratios.  It  is  also  evident  that  for  smaller  particles,  the  decrease  in  particle  size  shifts  the 
binodals  monotonically  toward  lower  concentrations  of  particles.  The  latter  is  indicative  of  more  extensive 
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immiscibility  for  smaller  particles  and  provides  strong  evidence  for  the  significance  of  the  multibody  effects 
incorporated  in  our  model 

Overall,  our  above  study  highlighted  several  important  aspects  of  depletion  in  nanoparticle-polymer 
mixtures  which  had  not  been  addressed  in  earlier  studies.  Explicitly,  considerations  of  the  depletion 
characteristics  (presented  in  ref  39)  suggested  that  the  asymptotic  results  derived  by  the  earlier  theories  are 
applicable  only  for  either  extremely  small  or  extremely  large  particles  and  that  crossover  effects  can  play  an 
important  role  in  determining  the  interactions  and  phase  behavior  of  intermediate- sized  particles.  More 
interestingly,  our  analysis  suggested  that  while  simple  considerations  based  on  surface  tension  effects  (Fig.  2) 
may  correlate  to  the  miscibility  at  dilute  concentrations,  the  overall  phase  behavior  and  dispersability  in  the 
nanoparticle  regime  may  exhibit  much  richer  characteristics.  In  the  next  section,  we  discuss  another  example 
where  similar  inferences  are  drawn. 

C.  Nanoparticles  with  adsorbing  polymers 

In  a  second  application  of  our  formalism,—  we  considered  the  phase  behavior  and  mechanical  properties  of 
nanoparticle-polymer  mixtures  for  which  the  particles  have  favorable  enthalpic  interactions  with  the  polymer. 
This  situation  constitutes  the  common  scenario  where  one  may  expect  to  achieve  “dispersion"  of  the  particles 
and  is  hence  of  significant  interest  for  polymer  nanocomposite  applications.  Yet  again,  our  research  was 
specifically  motivated  by  the  interplay  between  polymer-particle  interactions  and  the  particle  curvature  in 
determining  the  phase  behavior  and  the  particle  structure  in  such  situations. 

The  potential  richness  of  the  phase  behavior  and  the  impending  curvature  effects  becomes  manifest  when 
considering  both  the  polymer  concentration  and  particle  size  dependence  of  polymer- mediated  interparticle 
interactions  (at  infinite  dilution,  obtained  using  the  formalism  described  in  section  IIA)  displayed  in  Fig.  4a  and 
b.  For  the  case  of  a  particle  size  ratio  R/Rg  =  0.5,  we  observe  that  for  dilute  bulk  polymer  concentrations 
bulk?  the  interaction  potentials  are  monotonically  attractive  as  a  function  of  interparticle  distance.  At  such  dilute 
polymer  concentrations,  effects  arising  from  interpolymer  interactions  are  expected  to  be  relatively  weak. 
Consequently,  when  two  particles  are  brought  closer,  the  polymers  are  free  to  adsorb,  form  more  interparticle 
bridges  and  gain  energy  without  incurring  concomitant  entropic  costs.  These  effects  lead  to  a  strong, 
monotonic  attraction  between  particles  at  low  <Viik-  Upon  increasing  the  ambient  polymer  concentration,  the 
interactions  develop  a  nonmonotonic  character,  displaying  attraction  at  large  interparticle  separations  followed 
by  a  repulsive  behavior  at  smaller  interparticle  distances.  At  even  higher  bulk  concentrations,  it  is  seen  that  the 
interactions  become  monotonically  repulsive  with  the  interparticle  distance.  The  above  changes  in  the  character 
of  the  interparticle  potentials  may  be  rationalized  as  arising  from  an  increase  in  the  interpolymer  interactions 
arising  from  the  osmotic  confinement  effects  and  more  importantly  the  saturation  of  the  surfaces  by  the 
adsorbed  polymers.  Together  these  effects  lead  to  the  repulsive  interactions  noted  for  intermediate  and 
concentrated  situations. 


Fig.  4  (Figures  adapted  with  permission  from  ref  44).  (a)  Pair 
interaction  potentials  as  a  function  of  interparticle  distance  d 
(normalized  by  polymer  radius  of  gyration  Rg)  for  polymer-  to¬ 
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particle  size  ratio  R/Rg  =  0.5.  Bulk  concentrations  are  ^buik  = 

1.29  (■),  2.58  (♦),  3.87  (a),  5.16  (•)  where  Auik  represents 
the  bulk  polymer  concentration  normalized  by  the  overlap 
concentration.  The  inset  shows  corresponding  interparticle 
forces  as  a  function  of  interparticle  distance  d/Rg  for  <^,uik  = 

1.29,  2.58  and  3.87.  (b)  Corresponding  pair  interaction 
potentials  for  R/Rg  =  0.25. 

If  we  consider  the  influence  of  the  size  of  the  particle  upon  the  characteristics  of  the  pair  interaction 
potentials  ( cf.  Fig.  4b).  we  observe  that  the  interaction  potentials  become  weaker  with  decreasing  size  of  the 
particles.  To  rationalize  this,  we  recall  from  our  above  discussion  that  the  magnitude  of  the  attractive 
interactions  are  determined  by  the  number  of  bridging  segments  between  the  particles.  A  smaller  particle  has 
smaller  surface  area  leading  to  the  formation  of  lesser  number  of  bridges.  The  decrease  in  bridging  can  be 
argued  to  be  the  underlying  cause  of  the  overall  weaker  attraction  between  the  particles. 

It  is  interesting  to  note  that  the  above  results  contradict  the  conclusions  one  may  derive  based  on  the 
surface  tension  results  displayed  in  Fig.  2.  Indeed,  situations  of  favorable  polymer-surface  interactions  were 
expected  to  lead  to  negative  surface  tension  and  promote  dispersability.  In  contrast,  the  above  results  suggest 
that  such  situations  are  accompanied  by  strong  polymer-mediated  interparticle  attractions  and  potential 
immiscibility.  To  delineate  the  physics  of  such  polymer- nanoparticle  mixtures,  in  ref  43  we  translated  the 
above  interaction  potentials  to  phase  diagram  predictions  by  using  a  thermodynamic  theory  very  similar  to  the 
one  employed  in  the  context  of  depletion  interactions  (section  IIB).  Figs  5a  and  b  display  the  corresponding 
results  which  include  regions  of  immiscibility  as  delineated  by  the  fluid- fluid  coexistence  curves  for  particle 
sizes  R/Rg  =1.0  and  0.25.  For  all  the  particle  sizes,  we  observe  genetically  a  fluid  phase  at  low  polymer 
concentrations  (below  the  schematic  lower  boundaries  indicated  by  dashed  lines),  followed  by  bridging- 
induced  phase- separation  at  higher  polymer  concentrations  and  subsequently  a  stable  mixture  regime  at  even 
higher  polymer  concentrations.  The  latter  stabilization  arises  as  a  consequence  of  the  saturation  of  the 
adsorption  and  the  repulsive  interactions  at  the  higher  polymer  concentrations. 


(a)  (b) 

Fig.  5  (Figures  adapted  with  permission  from  ref  44).  Fluid-fluid  co-existence  curve  (open 
and  filled  symbols)  and  percolation  line  (solid  line  continued  as  a  dotted  line  into  the  region  of 
coexistence)  in  the  polymer  concentration  (</>overan)-particle  volume  fraction  ('/<)  plane.  The 
area  above  the  co-existence  curve  shows  the  one  phase  region  and  that  below  each  curve 
represents  the  two  phase  region  The  compositions  of  the  coexisting  phases  are  shown  by  tie 
lines  joining  the  “Hoc”  (filled  symbols)  and  supernatant  (open  symbols)  compositions.  The 
lower  boundary  for  the  two  phase  region  is  displayed  schematically  by  dash-dot  line.  Area  to 
the  left  ofpercolation  line  shows  fluid  phase  and  that  to  the  right  represents  percolated  phase: 
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(a)  R/Rg  =  1.0;  (b)  R/Rg  =  0.25. 

It  is  evident  by  comparing  Fig.  5a  and  b  that  the  relative  polymer-particle  sizes  play  an  important  role  in 
influencing  the  structure  and  phase  behaviors.  At  dilute  particle  concentrations,  we  observe  that  a  lowering  of 
the  R/Rg  ratio  shifts  the  upper  boundary  of  the  two  phase  region  to  lower  polymer  concentrations.  The  latter 
suggests  that  polymer-nanoparticle  mixtures  involving  smaller  particles  at  dilute  particle  concentrations  tend  to 
become  miscible  at  much  lower  polymer  concentrations  compared  to  the  larger  particles.  A  second  particle 
size  effect  is  observed  in  the  compositions  of  the  coexisting  phases  denoted  through  the  tie  lines  (shown  by  the 
lighter  lines).  We  observe  that  for  larger  particles,  the  concentrations  of  polymer  in  the  two  coexisting  phases 
are  more  or  less  the  same.  The  latter  suggests  that  the  phase  separation  in  such  systems  forms  a  particle-rich 
and  particle- depleted  phase,  both  rich  in  polymers.  On  the  other  hand,  for  smaller  particles,  the  phase 
separation  is  into  a  supernatant  phase  that  is  dilute  in  both  the  polymer  and  the  particles,  whereas  the  “floe" 
phase  is  rich  in  both  the  polymer  and  particles.  The  latter  trends  are  consistent  with  the  phenomena  of" 
complex  coacervation”  commonly  observed  in  the  context  of  protein-polysaccharide  mixtures.67"69 

In  many  applications  involving  nanoparticles,  such  as  in  their  use  as  rheological  modifiers,  more  detailed 
information  regarding  the  polymer-particle  floe  complexes  and  their  mechanical  properties  is  also  desired.67"69 
To  address  these  issues,  we  developed  a  new  simulation  strategy4 1’43’69,70  in  which  the  above  framework  was 
first  extended  to  deduce  the  probability  of  forming  polymer  bridges  between  two  particles.  Subsequently,  a 
semi-grand  canonical  Monte  Carlo  simulation  was  implemented  using  the  polymer- mediated  effective 
interparticle  potentials  while  mapping  out  the  connectivity  of  the  particles  using  such  bridging  probabilities. 
Knowledge  of  the  connectivity  allowed  us  to  quantify  the  cluster  sizes  of  polymer- bridged  particles  and  the 
resulting  mechanical  strength  of  such  complexes  through  simple  elastic  network  theories.  In  Fig.  6  we  display 
the  results  obtained  using  this  idea  which  provided  for  the  first  time  quantification  of  the  percolation, 
complexation  thresholds  and  mechanical  strengths  in  such  mixtures. 


(*)  (b) 

Fig.  6  (Figures  adapted  with  permission  from  ref  41):  (a) 

Percolation  probabilities  (P)  for  the  polymer- bridged  particle 
gels  displayed  as  a  function  of  particle  volume  fraction,  *1,  for 
different  particle- to-polymer  size  ratios,  R/Rg  =  2.0,  1.0,  0.5  for 
a  polymer  melt  matrix;  (b)  elastic  modulus  ( G ')  (deduced  using 
a  network  theory  for  elasticity)  as  a  function  of  particle  volume 
fraction,  *?,  for  particle- to-polymer  size  ratios,  RJRg  =  2.0,  1.0, 

0.5  for  a  polymer  melt  matrix. 

In  summary,  the  above- discussed  situation  again  provides  an  interesting  illustration  of  the  idea  that  single 
particle  miscibility  considerations,  while  useful,  may  not  necessarily  provide  a  complete  picture  of  the 
dispersability  and  the  resulting  structure  of  the  polymer- nanoparticle  mixture. 

D.  Applications  to  anisotropic  fillers 
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In  more  recent  studies  we  have  extended  our  formalism  to  study  the  effective  pair- interaction  potentials  and 
the  resulting  phase  behavior  and  percolation  transitions  of  nanorods  dispersed  in  solutions  of  adsorbing 
polymers.—  We  again  used  a  polymer  self-consistent  field  theory  (SCFT)  framework  in  conjunction  with  a 
Deijaguin  approximation  to  compute  the  polymer- mediated  orientation- dependent  pair  interaction  potentials 
between  cylindrical  nanorods.  A  modified  Flory  theory  and  a  simple  analytical  model  were  then  used  to 
deduce  the  different  equilibrium  phases  and  the  onset  of  percolation  for  nanorods  in  polymer  solutions.Lk77 
Yet  again,  rich  phase  behavioral  characteristics,  including  the  possibility  of  isotropic  and  nematic  phases,  were 
deduced.  We  delineated  results  quantifying  the  influence  of  polymer- surface  affinity,  polymer  concentrations, 
radius  of  rods,  and  aspect  ratio  of  rods,  on  the  topology  of  such  equilibrium  phases  and  percolation  regimes. 

E.  Other  related  theoretical  studies 

Considering  the  significant  practical  ramifications,  theoretical  models  and  computer  simulations  for  predicting 
the  phase  behavior  of  polymer-nanoparticle  mixtures  has  attracted  significant  interest.  Of  this  research,  the 
issue  of  “depletion”  interactions  in  polymer  solutions  has  specifically  involved  many  modeling  and  simulation 
studies  in  both  the  physics  and  the  materials  sciences  communities.  The  origin  and  nature  of  depletion 
interactions  were  first  elucidated  over  50  years  ago  by  Asakura  and  Oosawa  (AO),—  and  independently  by 
Vrij.—  The  AO  model  is  known  to  be  adequate  only  for  the  case  of  dilute  and  noninteracting  polymer 
solutions  and  only  when  the  radius  of  gyration  Rg  is  much  smaller  than  the  size  R  of  the  particle.  Subsequently, 
many  studies  have  examined  the  depletion  characteristics  in  particle-polymer  mixtures  for  situations  beyond 
these  limiting  constraints.  For  instance,  Meijer  and  Frenkel  used  a  lattice  simulation  approach  to  analyze  the 
case  of  dilute,  ideal  polymer  solutions  with  larger  polymer  to  particle  size  ratios,  and  showed  that  multibody 
interactions  between  the  particles  can  lead  to  significant  corrections  to  the  phase  behavior  predicted  by  the 
AO  theory.  Schweizer  and  coworkers76"81  pointed  out  the  breakdown  of  AO  theory  for  the  regime  Rg/R  ~ 
0(1),  and  have  subsequently  developed  integral  equation  approaches  which  incorporate  the  interactions 
between  the  polymers  and  treat  a  wide  range  of  polymer  and  particle  sizes.  Bolhuis,  Louis  and  co  workers 
have  developed  a  novel  coarse- graining  technique  to  treat  the  case  of  interacting  polymers  for  sizes  R/Rg  upto 
0(1  goth  Schweizer's  and  Bolhuis  et  al.'  s  researches  established  the  importance  of  interpolymer 
interactions,  and  also  delineated  the  resulting  phase  behavior  for  a  variety  of  parametric  conditions. 
Lekkerkerker,  Tuinier  and  coworkers66,84"87  have  pioneered  the  Gibbs  adsorption  and  free  volume  theories 
to  develop  the  phase  diagrams  for  mixtures  of  colloids  of  different  geometrical  shapes  in  both  ideal  and 
interacting  polymer  solutions.  Also,  alternative  approaches  invoking  perturbation  theories,—  and  cell  models— 
have  been  used  to  predict  the  phase  behavior  of  nanoparticle-polymer  solutions  characterized  by  depletion 
interactions. 

In  contrast  to  the  above  depletion  situation,  there  have  been  far  fewer  theoretical  studies  quantifying  the 
influence  of  favorable  enthalpic  interactions  between  the  polymer  and  the  particle  surface.  The  earliest  studies 
in  this  regard  were  based  on  scaling  theories  designed  to  address  the  adsorption  and  interaction  characteristics 
of  polymer  solutions  on  flat  plates  and  large  colloidal  particles.59,58,90,91  Various  molecular  modeling 
techniques  such  as  density  functional  theories^  and  integral  equation  theories—  were  also  successfully  applied 
to  extract  the  detailed  structural  descriptions  and  thermodynamic  properties  of  polymer  solutions  in  the 
presence  of  flat  surfaces.  Extensions  to  the  nanoparticle  regime  have  mainly  been  at  a  scaling  level  and  relate 
only  to  the  adsorption  characteristics  on  a  single  nanoparticle.^Mi 

More  pertinent  to  the  PNC  applications  are  studies  considering  the  interaction  and  phase  behavioral 
characteristics  of  nanoparticles  in  polymer  melts.—  The  earliest  studies  specific  to  the  polymer  melt  context 
came  from  the  groups  of  Gianellis—  and  Balazs  and  co  workers,—  who  used  polymer  self-consistent  field 
theory  to  analyze  the  thermodynamics  of  mixing  between  polymer  and  clay- like  fillers.  Their  results  pointed  to 
the  interplay  between  entropic  polymer  conformational  effects  and  polymer- filler  energetic  interactions  in 
controlling  the  equilibrium  state  of  polymer-  clay  mixtures.  Subsequently,  in  a  series  of  articles, -fL22.  Balazs  and 
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co workers  extended  the  above  studies  by  combining  polymer  self-consistent  field  theory  with  density 
functional  theories  to  predict  complete  phase  diagrams  for  mixtures  of  polymers  and  plate  and  rodlike  particles 
in  a  polymer  melt,  which  among  other  predictions  also  delineated  optimal  conditions  for  creating  stable 
dispersed  composites. 

Use  of  integral  equation  theories  for  addressing  the  above  issues  in  the  context  of  polymer  melts  has  been 
pioneered  by  Schweizer  and  coworkers.36’53’100’101  Such  theories  can  accommodate  a  finer  representation  of 
the  polymer  chains  relative  to  coarse-grained  elastic  thread  models  and  can  hence  yield  information  regarding 
polymer  packing  and  particle  structure  while  incorporating  the  multibody  particle  effects  (albeit,  in  an 
approximate  manner).  Specifically  Schweizer  and  coworkers  have  adapted  the  PRISM  theory101’53  to 
investigate  the  equilibrium  miscibility,  particle  dispersion  and  phase  separation  of  polymer-particle  mixtures  and 
polymer  nanocomposites.  They  have  used  this  formalism  to  examine  the  influence  of  particle  size,  degree  of 
polymerization  and  melt  density  upon  the  structure,  effective  forces  and  thermodynamics  of  polymer 
nanocomposites.  A  notable  prediction  from  such  theories  is  the  possibility  of  two  distinct  kinds  of  phase 
separation  behaviors  in  mixtures  of  particles  and  polymer  melts.  The  first  occurs  at  lower  monomer- filler 
attraction  strength  which  corresponds  to  an  entropic  depletion  attraction- induced  phase  separation.  The 
second  regime  occurs  at  a  higher  monomer- filler  adsorption  energy  and  involves  the  formation  of  an 
equilibrium  physical  network  phase  with  local  bridging  of  fillers  by  polymers.  Selected  applications  of  the 
PRISM  framework  to  PNCs  have  been  summarized  in  a  recent  review  article.— 

Direct  computer  simulation  approaches  have  also  been  used  to  address  the  equilibrium  structure  and  phase 
behavior  of  mixtures  of  spherical  nanofiller  units  dispersed  in  homopolymers.33’102"1 10  Specifically,  Smith  and 
coworkers—  and  Keblinski  and  coworkers—  used  molecular  dynamics  simulations  to  study  the  effective 
interactions  between  two  spherical  nanoparticles  in  a  polymer  melt.  These  studies  quantified  the  magnitude  of 
effective  interactions  as  a  function  of  different  physical  parameters  including  the  interaction  strength  between 
the  polymer  monomers  and  the  particle.  De  Pablo  and  coworkers 1 1 1  ■ 1 1 2  have  used  molecular  simulations  to 
analyze  the  depletion  interaction  and  deduce  the  effective  polymer- mediated  interactions  between  colloidal 
particles.  Their  results  were  consistent  with  earlier  theoretical  predictions  except  insofar  as  monomer  level  and 
chain  length  effects  not  captured  in  coarse-grained  mean- field  theories. 

III.  Mixtures  of  homopolymers  with  grafted  nanoparticles 

In  this  section,  we  briefly  review  some  of  the  developments  in  the  context  of  coarse-grained  modeling  of 
functionalized  nanoparticles  dispersed  in  homopolymer  melts.  As  noted  in  the  introduction,  functionalization 
of  particles  is  a  common,  and  in  many  cases  the  only  viable  strategy  to  facilitate  the  dispersion  of  nanoparticles 
in  polymer  matrices.  Experimentalists  have  explored  the  use  of  two  broad  classes  of  functionalizers  in  these 
efforts:  (i)  anionic  or  cationic  oligomeric  surfactants;26,30’113  and  (ii)  grafted  polymers.1 14-1 17  In  modeling  the 
effect  of  small  molecule  surfactants,  we  note  that  at  the  level  of  coarse-grained  modeling,  such  surfactants  are 
typically  at  the  same  scale  as  that  of  a  single  coarse-grained  polymer  unit.  Consequently,  as  a  first 
approximation,  the  conformational  degrees  of  freedom  of  the  surfactant  maybe  ignored,  and  phase  behavior 
of  nanoparticles  functionalized  by  such  surfactants  may  be  mapped  onto  the  behavior  of  polymer-particle 
mixtures  in  which  the  surfactant  effects  manifest  just  as  effective  enthalpic  interactions  between  the  polymer 
and  particle  (the  results  for  the  latter  situation  was  discussed  in  section  II C).  On  the  other  hand,  for  modeling 
situations  where  the  nanoparticles  are  grafted  with  either  a  polymer  or  longer  surfactants,  the  conformational 
degrees  of  freedom  of  the  grafted  functionalizer  need  to  be  accounted  for  in  quantifying  the  interactions  and 
phase  behavior  of  the  resulting  mixture.  In  the  following,  we  briefly  review  the  developments  that  have 
accompanied  the  latter  context.  Without  loss  of  generality,  we  refer  to  the  functionalizers  as  polymers  with  the 
understanding  that  they  may  equally  well  be  a  longer  molecule  surfactant. 

The  simplest  model  system  which  has  attracted  the  most  attention  from  both  an  experimental  and 
theoretical  perspective  is  one  where  the  grafted  polymer  is  chemically  identical  to  the  matrix  polymer  (termed 
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the  “autophobic”  case).  In  such  cases  there  are  no  competing  enthalpic  interactions,  and  the  polymer- 
mediated  interactions  and  phase  behavioral  characteristics  arise  primarily  from  the  entropic  effects  pertaining 
to  the  grafted  polymer  and  the  matrix  chains.  Much  of  the  theoretical  developments  in  this  area  have  grown  out 
of  the  seminal  descriptions  of  the  wetting  and  dewetting  of  polymer  melts  on  polymer  brushes  advanced  by 
Leihler  and  coworkers.118’119  Explicitly,  using  scaling  ideas  and  strong  segregation  theory  calculations  they 
delineated  the  regions  where  the  matrix  polymer  wets  or  dewets  the  brush.  In  the  former  case,  the  matrix 
chains  completely  penetrate  the  brush,  whereas  in  the  latter  case  the  matrix  chains  are  either  expelled  or 
penetrate  only  a  finite  zone  into  the  grafting  layer  (cf  Fig.  7).—  From  the  diagram  of  states  displayed  in  Fig. 
7c,  it  can  be  seen  that  the  overall  behavior  depends  on  the  degree  of  polymerization  of  the  grafted  chains,  N, 
the  chain  grafting  density  a,  and  the  degree  of  polymerization  of  the  free  host  chains,  P.  Explicitly,  when 

try  \  -H  \  l  )  dewetting  0f  the  melt  chains  is  expected  (“dry”  brush).  In  other  words,  situations  involving 
long  matrix  polymers  and/or  densely  grafted  polymer  functionalizers  would  be  more  conducive  to  dewetting. 
We  also  note  that  while  these  considerations  were  derived  based  on  scaling  arguments  and  analytical  theories, 
Matsen  and  Gardiner  have  carried  out  careful  SCFT  based  numerical  analysis  of  these  issues.—  Their  results 
have  confirmed  many  of  the  above  details  at  a  qualitative  level  with  however  some  significant  quantitative 
differences  which  highlight  the  approximations  inherent  in  analytical  theories. 


(a) 


(b) 
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Fig.  7  (a)  and  (b)  A  schematic  of  the  wetting  to  de  wetting  transition  in  the  interlaces  between  melt  and  polymer 
brushes.  The  wetting  regime  corresponds  to  the  penetration  of  the  melt  chains  into  the  brush.  In  contrast,  the 
dewetting  regime  corresponds  to  either  the  partial  penetration  or  complete  expulsion  of  the  melt  chains  from  the 
brush.  These  behaviors  are  indicated  schematically  also  in  the  volume  fraction  profiles  of  the  melt  (M)  and  brush 
chains  (B).  (c)  (Reproduced  with  permission  from  ref  119)  Diagram  of  wetting-dewetting  transitions  in  ( N ,  P,  cf) 
plane  for  a  polymer  brush  of  polymerization  index  A  immersed  in  a  melt  of  the  same  polymer  with  a  different 
polymerization  index  P.  The  full  lines  are  the  boundaries  between  the  regions  with  different  scaling  laws  for  the  brush 
height,  h.  The  dashed  line  separates  the  two  regimes  of  scaling  of  the  brush- melt  interfacial  thickness  and  coincides 
with  the  frontier  between  the  regions  of  positive  and  negative  spreading  coefficients. 
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To  translate  the  above  considerations  to  the  dispersion  of  nanoparticles,—  we  observe  that  the  case  of  “ 
dewetting”  is  akin  to  the  polymer  “depletion”  case,  and  hence  one  may  expect  that  for  extreme  cases  of 
dewetting  the  matrix  polymer- mediated  interparticle  interactions  become  attractive,  and  lead  to  aggregation 
and  phase  separation  of  the  particles.  In  contrast,  the  “wetting”  situation  is  similar  to  one  where  the  matrix 
polymers  possess  favorable  enthalpic  interactions  with  the  particles  (except  insofar  as  the  absence  of  polymer 
bridges),  and  hence  the  interparticle  interactions  and  phase  behavior  may  be  expected  to  favor  dispersion  or 
mixing.  Such  considerations  were  first  confirmed  in  a  combined  theoretical  and  experimental  work  by 
Hasegawa  and  coworkers  who  calculated  the  interplay  between  the  brush-brush  repulsions  and  the 
emergence  of  dewetting- induced  attractive  interactions.—  An  interesting  optimum  intermediate  grafting  density 
was  predicted  where  the  net  attractions  were  weakest  and  the  particles  are  most  well-dispersed.  More 
recently,  systematic  experiments  by  Green  and  coworkers  and  others  have  confirmed  the  correspondence 
between  the  wetting-dewetting  transitions  and  the  miscibility  behavior  of  grafted  nanoparticles. 1 22~ 1 24 

While  much  of  the  above- discussed  studies  were  based  on  results  obtained  by  considering  the  behavior  of 
flat  grafted  surfaces  (in  conjunction  with  possibly  Deijaguin  approximations),  very  recently  density  functional 
theories,125’126  SCFT  approaches127’128  and  computer  simulations129"132  have  started  tackling  the  particle 
curvature  effects  in  a  much  more  direct  manner.  For  instance,  the  interaction  between  two  brush  coated 
spheres  in  a  good  solvent  has  been  studied  to  confirm  its  purely  repulsive  nature.128’125  Other  density 
functional  studies  and  simulations125’133  have  also  examined  the  influence  of  solvent  quality  effects,  and  have 
suggested  that  richer  interaction  characteristics  including  attractions  and  repulsions  are  possible  depending  on 
the  solvent  quality.  Recent  SCFT  studies  for  grafted  nanoparticles  in  polymer  melts  have  also  examined  in 
more  detail  the  effects  of  particle  size  and  grafting  density  upon  the  interparticle  interactions  (albeit,  for  a 
specific  ratio  of  grafted  and  free  polymers  chain  lengths).—  The  latter  studies  point  to  the  applicability  of 
considerations  gleaned  from  the  studies  on  flat  plates,  with  however  much  weaker  interactions  and  a  particle 
size,  grafting  density  dependent  shift  of  the  wetting-dewetting  regimes  for  curved  surfaces.  The  implications  of 
such  weaker  interactions  for  the  context  of  dispersion  of  nanoparticles  have  been  examined  using  a  simple 
theory  by  Harton  and  Kumar.— 

Another  independent  line  of  investigation  has  been  initiated  by  Schweizer  and  coworkers  focusing  on  the 
use  of  integral  equation  theories  to  shed  light  on  the  behavior  of  “sparsely  grafted”  nanoparticles  (with  only  a 
few  grafted  chains)  dissolved  in  a  homopolymer  matrix.36’136"139  PRISM  theories  have  gone  beyond  many 
above-mentioned  studies  by  considering  the  effects  arising  from  finite  concentration  of  particles.  Interestingly, 
they  predict  that  melts  and  dense  solutions  of  nanoparticles  may  exhibit  signatures  of  “microphase  ordering” 

(in  addition  to  the  possibility  of  macrophase  demixing)  accompanied  by  the  formation  of  structural 
characteristics  at  a  finite  length  scale. 

A.  Some  future  directions 

In  concluding  this  section,  we  mention  two  issues  which  have  received  far  less  attention  in  the  context  of 
dispersing  polymer  functionalized  nanoparticles  in  polymer  melts: 

(i)  A  complete  theoretical  understanding  of  the  interactions  and  phase  behavior  of  polymer  grafted 
nanoparticles  dispersed  in  homopolymer  melts  requires  consideration  of  a  vast  parameter  space  in  which  the 
molecular  weight  of  the  grafted  polymer  and  its  chemical  identity  need  to  be  accounted  in  addition  to  effects 
arising  from  the  particle  size,  the  particle's  interactions  with  the  polymer  matrix  and  the  molecular  weight  and/or 
concentration  of  the  matrix  polymer.  The  studies  discussed  in  the  preceding  section  have  clarified  the  roles  of 
the  relative  sizes  of  the  matrix  and  grafted  polymer  lengths,  grafting  densities  and  (to  a  more  limited  extent)  the 
curvature  of  the  particles.  However,  an  issue  which  has  attracted  less  attention  has  been  the  role  of  chemical 
mismatch  between  the  matrix  and  grafted  polymers.  Indeed,  the  parameter  space  available  to  synthetic 
chemists  is  considerably  enhanced  if  polymers  chemically  distinct  from  the  matrix  polymer  are  used  to 
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functionalize  and  disperse  the  particles. 

A  seminal  study  in  the  above  regard  was  carried  out  by  Borukhov  and  Leibler,  where  the  tethered  and 
matrix  polymers  were  allowed  to  be  chemically  distinct  but  with  favorable  interactions  (the  brush- matrix  Flory 
interaction  parameter  was  chosen  to  be  negative).140’141  They  used  scaling  arguments  to  show  that  under 
appropriate  conditions  this  may  serve  to  eliminate  the  effective  attraction  (and  potential  immiscibility)  noted  in 
the  case  where  the  matrix  polymers  are  chemically  identical  to  the  grafted  chains.  More  recently,  we  have 
examined  the  wetting  characteristics  of  polymers  which  are  chemically  different  and  possess  unfavorable 
interactions  with  the  brush  component.—  We  reported  experimental  results  and  complementary  strong- 
segregation  theory  arguments  on  the  parametric  interplay  between  enthalpic  and  entropic  effects  in  the 
interfaces  between  polymers  and  polymer  brushes.  Our  studies  indicated  (cf.  Fig.  8)  that  one  may  be  able  to 
use  brushes  made  of  incompatible  polymers  of  sufficiently  high  molecular  weight  and  achieve  lower  interfacial 
tensions  compared  to  brushes  made  of  compatible  polymers  of  small  molecular  weight.  Since  overall  particle 
dispersability  usually  correlates  to  the  melt- brush  interfacial  tensions,  this  strategy  may  open  the  door  to  more 
functionalization  possibilities  when  synthesis  and/or  grafting  methods  prove  to  be  limiting. 


Brush  thickness  by  PSOH  (inn) 

Fig.  8  (Reproduced  with  permission  from  ref  142) 
Strong  segregation  theory  calculations  of  interiacial 
tension  I'ps-pmma  between  polystyrene  brush  and  a 
poly(methyl  methacrylate)  (PMMA)  melt  for  different 
Mn  of  the  PMMA  polymer  melt.  Po  and  b 
respectively  denote  the  monomeric  volume  and  the 
segment  length  of  the  polymer.  The  closed  symbols 
were  computed  using  a  Flory-Huggins  interaction 
parameter  X  =  /ps-pmma  =  0.037.  The  open 
symbols  used  X  =  0  which  corresponds  to  the  “ 
autophobic”  case.  It  is  evident  that  in  practice  one 
maybe  able  to  use  brushes  made  of  incompatible 
polymers  of  sufficiently  high  molecular  weight  and 
achieve  lower  interiacial  tensions  compared  to 
brushes  made  of  compatible  polymers  of  small 
molecular  weight. 


(ii)  A  second  issue  relates  more  closely  to  the  theme  of  the  studies  discussed  in  the  previous  sections,  viz., 
behavior  arising  at  nondilute  concentrations  ofnanoparticles.  Indeed,  most  of  the  studies  mentioned  above 
(except  the  recent  PRISM  efforts)  relate  to  either  the  wetting/dewetting  considerations  or  the  interactions 
arising  in  the  context  of  two  particles.  While  such  results  provide  valuable  guidelines  for  dispersion  strategies, 
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the  structure  resulting  in  the  multiparticle  situation  may  potentially  exhibit  much  richer  features.  A  recent 
example  of  this  effect  was  noted  in  the  context  of  experiments  and  related  theoretical  studies  on  dispersing 
spherical  nanoparticles  grafted  with  polymeric  brushes  into  a  homopolymer  matrix  (see  Fig.  9). 143444  The 
theoretical  ideas  were  based  on  simple  scaling  models  to  enumerate  the  free  energies  of  the  different  structures 
observed  in  computer  simulations  (and  experiments).  Explicitly,  it  was  proposed  that  the  equilibrium  structures 
were  chosen  as  a  result  of  the  competition  between  the  entropic  penalty  arising  due  to  the  conformational 
rearrangements  of  the  grafted  chains  and  the  enthalpic  gain  arising  from  the  bare  particle-particle  attractions. 
Even  within  the  “wetting”  regime,  we  noted  that  the  preceding  interplay  between  can  lead  to  novel  self- 
assembly  of  the  particles  into  anisotropic  structures.  While  this  study  was  just  an  isolated  example  highlighting 
the  issues,  the  rich  characteristics  of  the  multiparticle  assembly  and  phase  behavior  for  the  full  parameter  space 
of  dispersion  of  polymer- grafted  particles  in  polymer  melts  still  remains  to  be  elucidated. 


0  4  8  12  IS 

Grafted  Child  length 


G  rifled  chi  in  lerglh  t  matrix  chain  length 


Fig.  9  (Adapted  with  permission  from  ref  143)  Parametric 
phase  diagram  for  the  structures  formed  during  the  dispersion  of 
polymer  grafted  nanoparticles  in  a  polymer  matrix.  The  matrix 
polymer  was  kept  the  same  while  the  number  and  size  of  the 
grafted  polymers  were  varied,  (a)  A  comparison  of  theoretic 
calculations  based  on  strong- segregation  theory  (solid  lines)  and 
simulations  (points).  Spherical  symbols:  spheres,  square 
symbols:  sheets,  triangles:  strings,  diamonds:  well- dispersed 
particles;  (b)  Experimental  ‘morphology  diagram'  of  polymer- 
tethered  particles  mixed  with  matrix  polymers.  Spherical 
symbols:  spheres,  square  symbols:  sheets,  triangles:  strings, 
diamonds:  well- dispersed  particles.  The  lines  that  separate  the 
different  regions  are  merely  guides  to  the  eye  (see  ref  143  for 
details  on  the  experimental  system). 


IV.  Self-assembly  of  nanoparticles  in  block  copolymer  matrices 

Recently,  the  lessons  learned  from  efforts  to  maximize  dispersion  of  nanoparticles  in  homopolymer  melts  have 
been  furthered  towards  “dispersion  control”  which  focuses  on  directing  the  self-assembly  in  nanoparticle- 
block  copolymer  mixtures.  These  ideas  have  led  to  a  surge  of  experimental  reports  which  exploit  the  self- 
assembly  in  mixtures  of  diblock  copolymers  and  nanosized  particles  to  produce  ordered  organic- inorganic 
hybrid  materials.7’145"151  In  some  applications,  the  microphase  separation  of  the  block  copolymers  is  used  as 
a  template  to  control  the  ordering  of  the  particles  and  to  produce  highly  organized  hybrid  materials.145"147’152 
The  possibility  of  using  such  strategies  to  achieve  significant  loading  of  nanoparticles  has  also  been 
demonstrated.—  In  other  applications,  the  particles  are  used  to  modify  the  self-assembly  of  the  parent  block 
copolymer  to  lead  to  new  morphologies  of  self-assembly.154’155  The  resulting  structures  have  been  proposed 
for  use  in  applications  such  as  separation  processes,  next- generation  catalysts  and  photonic  band  gap 
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materia  Is. Til? 

For  successful  fruition  of  the  above  applications,  a  fundamental  understanding  of  the  manner  in  which 
different  parameters  in  such  systems,  such  as  the  size,  shape,  volume  fraction  of  particles,  copolymer 
composition,  and  interaction  energies  between  the  different  components  control  the  thermodynamics  and  self- 
assembly  of  such  nanoparticle-block  copolymer  mixtures.  Nanoparticle-block  copolymer  composites 
represent  another  facet  in  the  category  of  nanoparticle-polymer  mixtures  where  simple  concepts  such  as  the 
surface  tension  effects  ( cf.  Fig.  2  and  the  accompanying  discussion)  alone  are  not  expected  to  provide  a 
complete  understanding  of  the  richness  of  the  phase  behavior  and  dispersion  characteristics.  Indeed,  the 
overall  self-assembly  in  such  systems  is  expected  to  depend  on  an  intricate  interplay  of  such  surface  tension 
effects  with  the  energetic  effects  driving  the  block  copolymers  to  self- assemble.  In  this  section,  we  briefly 
review  some  of  the  theoretical  developments  which  have  occurred  used  coarse-grained  modeling  and 
simulation  to  clarify  these  effects  in  the  context  of  nanoparticle  organization  and  self-assembly  in  block 
copolymers. 

A.  Templated  organization  of  nanoparticles  in  self-assembled  phases  of  block  copolymers 

In  many  applications,  it  is  desired  to  achieve  templated  organization,  in  which  the  nanoparticles  are  either 
directed  to  the  interlace  or  to  exist  wholly  within  one  of  the  phases  of  the  self- assembled  block  copolymer 
phases.  A  fundamental  question  confronting  such  strategies  is:  “what  are  the  physical  parameters  controlling 
the  nanoparticle  distributions  in  block  copolymers?”  This  question  was  first  addressed  using  modifications  of 
SCFT  theories  by  Balazs,—  and  then  subsequently  by  using  molecular  dynamics,—  Monte  Carlo 
framework!^-  and  a  hybrid  field  theory  based  simulation  approach.^  Broadly,  the  results  of  these  studies 
suggested  that  the  templating  of  the  particles  by  the  block- copolymer  is  dependent  on  the  size  of  the  particles 
and  their  interactions  with  the  different  units  of  the  copolymer.  If  the  particles  were  compatibilized  to  just  one 
of  the  components  (“selective”  particles),  then  they  were  predicted  to  localize  at  the  center  of  their  preferred 
phase,  while  particles  compatible  to  both  components  (“nonselective”  or  “surfactant- like”  particles)  were 
predicted  to  localize  at  the  AB  interlace  of  a  AB  diblock  copolymer. 

In  recent  work,  we  used  strong- segregation  approximation  to  develop  an  analytical  theory  which  provides 
a  mechanistic  basis  and  identifies  the  important  parameters  governing  the  above  results  of  particle 
distributions.—  Explicitly,  we  argued  that  there  were  three  primary  energetic  factors  whose  interplay  governed 
the  particle  distributions:  (i)  particles  positioned  at  the  interlace  of  the  copolymeric  phases  decrease  the 
interfacial  contacts  between  the  blocks  and  lower  the  accompanying  interfacial  energy  costs;  (ii)  particles 
positioned  in  their  preferred  phase  (i.e.  the  phase  to  which  the  particles  possess  relatively  more  favorable 
enthalpic  interactions),  gains  in  energy;  and  (iii)  location  of  the  particle  within  the  brush- like  block  copolymer 
phases  incurs  elastic  energy  costs  arising  from  the  distortion  of  the  brush  chains.  These  elastic  costs  are 
expected  to  be  largest  in  magnitude  near  the  interface  of  the  blocks  (the  “grafting”  location)  and  weakest  at 
the  locations  which  are  furthest  from  the  interlaces. 

To  render  the  above  arguments  concrete,  in  the  following,  we  consider  the  lamellar  phase  of  a  symmetric 
AB  diblock  copolymer  and  denote  the  particle  radius  as  R.  We  take  a  coarse-grained  view  where  the 
interactions  between  the  particle  and  the  A  and  B  components  of  the  polymer  are  quantified  respectively  by 
two  interfacial  tension  parameters  denoted  *]Ac  and  *lBc-  The  parameter  6*1  =  (*1ac  ~  ’ he )  represents  the  “ 
selectivity”  of  the  particle  to  the  polymer  component  (for  S'}  <  0,  the  particles  are  preferential  to  the  A  phase). 
By  denoting  >}AB  to  be  the  interfacial  tension  between  the  A  and  B  phases,  the  energetic  terms  (i)  and  (ii) 
above  can  be  estimated  as  a  function  of  the  nanoparticle  location  z  ( cf.  Fig.  1 0a)  relative  to  the  AB  interfacial 
plane  as 
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Fig.  10  (a)  A  schematic  of  the  particle  configuration  relative  to  the  AB  interlace  (denoted  by  the 
plane  z  =  0).  (Adapted  with  permission  from  ref  42):  (b)  and  (c)  Probability  distribution  of  the 
nanoparticle  location  in  the  block  copolymer  lamella  as  a  function  of  the  selectivity  parameter  ft 
and  the  particle  size  R  in  Rg  units.  In  (b)  the  particle  size  R  =  0.5.  z  denotes  the  distance  normal  to 
the  plane  of  the  AB  interlace  (in  Rg  units)  with  z  =  0  corresponding  to  the  AB  interlace.  z= *  1 
corresponds  to  the  center  of  the  A  layer). 


In  the  above,  the  two  distinct  cases  arise  from  the  reduction  in  the  AB  interlace  when  the  particle  is 
positioned  such  that  \z\  <  R,  and  the  lack  of  such  an  effect  for  \z\  >  R.  The  term  (iii)  above  can  be  estimated 
using  the  results  of  Williams  and  Pine  us—  and  the  more  recent  work  of  O'Shaughnessy  and  Kinf-^  as 


A E*tz)=P{z) 


3 


(7) 


where  P(z)  denotes  the  osmotic  pressure  field  acting  on  the  segments  of  a  polymer  brush  at  a  location  z.  Using 
eqn  (6)  and  (7),  we  can  approximate  the  density  distribution  of  nanoparticles  Pcftz)  in  the  diblock  lamella  as: 

Pc(z)  CC  exp  [— F enth(z)  -  AFbr(z)].  (8) 


By  using  expressions  from  strong  segregation  theory  for  P(z)  and  *Iab  we  can  deduce  the  parametric 
dependencies  of  the  density  distribution  of  the  nanoparticles.42-159  In  the  following,  we  briefly  discuss  ( cf.  Fig. 
10b  and  c)  the  results  obtained  as  a  function  of  the  parameter  ft  =  3f1/flAB  (ft  1  quantifies  the  degree  of 
amphiphilicity  or  the  surfactant- like  nature  of  the  particle)  and  the  size  R.  From  Fig.  10b.  it  can  be  seen  that  for 
a  fixed  R,  small  values  of  ft  lead  to  a  localization  of  the  particles  at  the  AB  interlace  of  the  copolymer.  In 
contrast,  increasing  the  selectivity  of  the  particles  leads  to  a  delocalization  of  the  particles  into  their  preferred 
phase.  Interestingly,  we  observe  that  for  intermediate  values  of  selectivity,  the  overall  density  distribution 
displays  three  peaks  corresponding  to  a  localization  at  the  middle  of  the  brush  in  its  preferred  phase.  The  latter 
is  a  manifestation  of  the  interplay  between  the  lower  elastic  energetic  cost  associated  with  the  particle  being 
present  at  the  top  extremities  of  a  polymer  brush  compared  to  its  interiors  and  the  interfacial  energy  gain  in 
localizing  at  the  interface. 

Size  effects  are  presented  in  Fig.  10(c).  It  can  be  seen  that  small  selective  particles  are  predicted  to  be 
more  localized  at  the  AB  interlace,  whereas  the  larger  particles  tend  to  exhibit  more  preferential  segregation. 
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These  effects  can  be  rationalized  as  arising  from  the  feet  that  the  interfacial  tension  gain  (eqn  (6))  scales  as  R2 
(surface  area)  whereas  the  elastic  energy  cost  (eqn  (7))  scale  as  R3  (the  volume  of  the  particle).  Consequently, 
for  small  particles  and/or  for  stronger  segregation  between  A  and  B  phases  (i.e.  larger  AB  interfacial  tension), 
the  particles  can  be  expected  to  be  more  localized  at  the  AB  interface.  In  contrast,  for  larger  particles  and/or 
weaker  segregations,  the  tendency  to  segregate  into  the  preferred  domain  dominates. 

A  related  outcome  of  the  above  analysis  was  the  prediction  that  the  addition  of  “surfactant- like" 
nanoparticles  (i.e.  with  selectivity  Stf  0)  are  expected  to  contract  the  lamellae  and  lower  the  elastic 
constants  of  the  block  copolymer.  Both  these  effects  can  be  physically  understood  as  arising  from  the 
reduction  in  the  AB  interfacial  costs  arising  from  the  positioning  of  the  particles  at  the  interface.  Hence,  the 
chains  have  to  stretch  less  to  accommodate  the  unfavorable  AB  contacts.  A  quantitative  analysis  of  such 
effects  also  suggested  that  lowering  of  the  elastic  modulii  may  lead  to  the  nanoparticle- induced  creation  of 
bicontinuous  phases  in  the  block  copolymer.  Shown  in  Fig.  1 1  are  experimental  results  confirming  such 
predictions.— 


<a)  <b) 

Fig.  11  (Adapted  with  permission  from  ref  154)  (a)  The  volume  fraction  (^)  dependence  of 
lamellar  thickness  h($)  (normalized  by  the  values  at  zero  particle  concentration  h0)  for  PS-6- 
P2VP  diblock  copolymers  with  M„  values  of  1 14  kg  mol-1  (squares),  196  kg  moF1  (circles), 
and  380  kg  mol-1  (triangles).  The  lines  (M„  values  indicated)  correspond  to  the  predictions  of 
strong  segregation  theory,  (b)  Cross-sectional  TEM  images  ofPS-/?-P2VP  block  copolymer 
(M„  196  kg  mol-1)  containing  PS-Au  nanoparticles  at  a  volume  fractions  of  0.09.  It  is  evident 
that  the  microstructure  of  the  PS  and  P2VP  domains  becomes  bicontinuous. 


In  closing,  we  note  that  the  above  considerations  were  based  on  analytical  arguments  founded  on  strong- 
segregation  theory  calculations.  More  recently,  Kim  and  Matsen  have  presented  a  careful  quantitative  analysis 
of  the  particle  distributions  using  a  novel  numerical  implementation  of  the  SCFT  formalism160"162  While  then- 
results  for  the  bare  particles  are  qualitatively  consistent  with  the  arguments  presented  in  ref  42,  they  have  also 
extended  these  considerations  further  by  treating  accurately  the  influence  of  polymeric  functionalizers. 

B.  Self-assembly  in  block  copolymer  nanoparticle  composites 

The  developments  discussed  in  the  preceding  section  pertain  to  the  physics  of  templated  assembly  of 
nanoparticles  in  block  copolymer  matrices.  However,  a  frill  understanding  of  the  morphology  of  the  block 
copolymer- nanoparticle  composites  must  accommodate  the  possibility  of  both  particle  self-assembly  as  well 
as  particle- induced  modifications  of  the  block  copolymer  self-assembly.  Seminal  steps  towards  a  theoretical 
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description  of  this  problem  was  taken  by  Balazs  and  coworkers,  who  extended  the  self-consistent  field  theory 
of  multicomponent  polymers  to  include  the  presence  of  hard  particles  of  different  shapes  and  delineated  the 
resulting  particle  and  block  copolymer  self-assemblies  for  a  variety  of  physical  parameters  which  included 
confinement  effects.37’ 163-1 67  Broadly,  the  results  of  their  analyses  suggested  a  rich  self-assembly  behavior 
determined  by  an  interplay  between  the  shape,  size  and  selectivity  of  the  particles  and  the  other 
physicochemical  features  of  the  block  copolymers.  Recently,  molecular  dynamics  simulations,—  cell  dynamics 
based  approaches,—  density  functional  theories—  and  Monte  Carlo  simulations170,171  have  also  been  used 
to  study  similar  issues.  These  studies  have  suggested  phase  behavior  that  is  qualitatively  consistent  with  the 
predictions  of  the  SCFT  theory  of  Balazs  and  coworkers. 

C.  Some  future  directions 

While  a  number  of  advances  have  been  reported  recently  in  the  theoretical  modeling  and  simulations  of  the 
interplay  between  nanoparticles  distributions  and  the  self-assembly  in  block  copolymer  phases,  a  question 
which  is  yet  only  partially  resolved  is  the  impact  of  surfactants  and  polymeric  functionalizers  in  the  assembly  of 
nanoparticles  in  block  copolymer  phases.  In  situations  where  the  functionalizers  are  small  molecule  surfactants, 
their  influence  may  be  subsumed  within  effective  energetic  parameters,  and  models  reviewed  in  the  preceding 
sections  may  suffice  to  identify  the  parameters  controlling  the  particle  distribution  and  block  copolymer  self- 
assembly.  However,  the  more  interesting  and  practically  important  case  is  one  where  the  nanoparticles  contain 
grafted  polymers.  As  mentioned  above,  Kim  and  Matsen—  recently  presented  an  analysis  of  such  a  single 
particle  case  to  deduce  the  distribution  of  such  nanoparticles  in  the  block  copolymer  phases.  In  earlier  studies, 
Balazs  and  co  workers—  considered  the  case  of  nondilute  concentrations  of  particles  but  each  containing  just 
one  grafted  polymer  (referred  to  as  a  “tadpole"  configuration).  They  extended  their  density  functional  theories 
to  address  the  self-assembly  in  such  cases.  In  other  research,  Reister  and  Fredrickson!^-  used  a  creative  idea 
of  modeling  the  grafted  nanoparticle  as  a  star  polymer  with  a  finite  sized  (soft)  core  to  shed  light  on  the  self- 
assembly  behavior  one  might  expect.  While  these  preceding  studies  and  their  results  have  been  invaluable, 
issues  such  as  the  role  of  the  molecular  weight  of  the  grafted  polymer  (relative  to  the  matrix  molecular  weight), 
the  grafting  density  of  the  nanoparticles  and  enthalpic  interactions  (if  any)  between  the  matrix  polymer  and  the 
functionalizers  are  still  unresolved  and  are  expected  to  constitute  active  directions  for  future  theoretical 
research. 

V.  Concluding  remarks 

In  summary,  we  briefly  reviewed  some  of  the  recent  theoretical  developments  in  the  context  of  coarse-grained 
modeling  of  equilibrium  characteristics  of  particle  dispersion  in  homopolymer  and  block  copolymer  matrices. 
The  studies  reviewed  were  connected  thematically  by  pointing  out  that  in  each  of  the  instances,  knowledge  of 
particle-polymer  interactions  at  the  single  particle  level  may  not  alone  suffice  to  explain  the  dispersion  and 
organization  characteristics  of  the  nanoparticles.  In  each  case,  we  highlighted  our  contributions  to  the  specific 
problem  at  hand  and  mentioned  some  related  theoretical  research  and  future  directions.  Admittedly,  our 
perspective  was  biased,  not  only  in  the  topics  reviewed  but  also  in  the  emphasis  on  certain  class  of  modeling 
approaches. 

We  note  that  our  discussions  exclusively  focused  on  theoretical  descriptions  of  the  “equilibrium” 
characteristics  of  particle  dispersion  and  assembly.  However,  despite  the  best  experimental  strategies, 
nonequifbrium  effects  resulting  from  spin  casting,  filler  aggregation  and/or  external  fields  are  bound  to  remain 
important  for  many  applications  ofPNCs.1 16,174  In  this  regard,  many  unresolved  theoretical  questions  still 
remain:  “how  does  the  structure  of  a  PNC  dispersion  evolve  upon  dispersing  the  fillers  in  the  polymer  matrix?” 
,  “can  quantities  such  as  the  cluster  size  distributions  and  fractal  dimensions  be  predicted  for  specified 
polymer- filler  combinations?",  “how  does  externally  applied  shear,  electric  and  magnetic  fields  (and 
combinations  thereof)  impact  upon  the  nonequilibrium  state  of  the  dispersion?”32,175  While  traditional 
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computer  simulations  may  shed  light  on  some  of  the  relevant  issues  pertaining  to  these  questions,  the  time  and 
length  scales  which  can  be  probed  by  such  means  may  not  necessarily  overlap  with  experimental  regimes,  and 
there  is  a  need  for  development  of  new  approaches  to  address  the  pertinent  issues. 

We  emphasize  that  coarse-grained  modeling  is  but  one  rung  in  the  ladder  of  modeling  approaches  for 
materials  structure  and  properties.  Other  approaches  such  as  quantum  mechanical  calculations,  atomistic 
simulations,  integral  equation  theories  and  continuum  mechanical  approaches  provide  complementary 
information  to  effect  predictive  computer  modeling  of  the  structure  and  properties  of  materials.  This 
complementarity  becomes  most  evident  when  one  desires  to  relate  the  parameters  accompanying  coarse¬ 
grained  models  to  the  chemical  details  of  the  polymer  molecules  and  the  filler.  This  requires  the  development  of 
efficient  multiscale  computer  simulation  tools  and  methodologies  which  can  render  quantitatively  the  connection 
between  the  chemistry  of  the  components,  their  force  fields  and  the  coarse-grained 
parameters.34’35’49’51’176,177  Availability  of  such  a  suite  of  tools  will  render  the  ah  initio  computer-aided 
predictive  characterization  of  properties  ofPNCs  a  reality. 
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